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Identification of r leptons at the D0 experiment 

R. Madar a for the D0 Collaboration 
a CEA/DSM/IRFU/SPP Saclay, FRANCE 

The article describes the identification of hadronically decaying r leptons in pp collisions at 1.96 TeV collected 
by the D0 detector at the Fermilab Tevatron. After a brief description of the motivations and the challenges of 
considering r leptons in high energy hadronic collisions, details of the r reconstruction and identification will be 
discussed. The challenges associated for r energy measurements in an hadronic environment will be presented 
including approaches to deal with such measurements. 



1. MOTIVATIONS AND CHALLENGES 

Leptons are of particulary interest in high 
energy hadronic collisions because they provide 
clean event signatures in a complex hadronic en- 
vironment. Electrons and muons are often consid- 
ered because of such detector signatures. In addi- 
tion, r leptons (called simply r here after) brings 
a potential acceptance gain of 50 — 300 % accord- 
ing to the lepton multiplicity of the analyzed final 
state. Moreover, studies with r leptons provide 
an interesting area to test the consistency of the 
Standard Model (SM) through lepton universal- 
ity at high energies by measuring the branching 
ratio of electroweak (EW) bosons decaying into 
r leptons [1 . One promising candidate for ex- 
tension to the SM is the Minimal Supersymmet- 
ric Standard Model (MSSM) which predicts new 
particles that can decay into r leptons. The exis- 
tence of such scenario can be probed within r final 
states |2|3j . Finally, many decay chains initiated 
by the Higgs boson involve r leptons and con- 
sidering them can then increase the experimental 
sensitivity to understand the origin of electroweak 
symmetry breaking (EWSB) [4]. 

Considering r leptons thus appears to encom- 
pass many virtues. However, understanding such 
an object in a hadronic environment is very chal- 
lenging. The neutrino (s) from r final states es- 
cape the detector without interacting and forces 
a fraction of the r energy to be invisible. Hence, 
the resulting visible energy is relatively soft and 
includes contamination from soft QCD interac- 



tions, which becomes a primary component of 
the hadronic final state. Further, the various 
decay mode of r must be considered. On the 
one hand, the identification of r from its lep- 
tonic decay {BIZ ~ 35%) suffers from electrons 
and muons originating from direct EW bosons de- 
cay as well as a poor statistic due to B1Z{tt — >• 
e/i) = 6%. On the other hand, r identification 
from its hadronic decays {BIZ ~ 65%) have differ- 
ent detector signature according to the hadronic 
resonance involved. Moreover, direct QCD inter- 
actions from hadrons collisions produce a lot of 
hadronic final states which can mimic the hadron- 
ically decaying r. To be optimal, each of the dif- 
ferent hadronic r final states requires dedicated 
analysis and a combination is performed to ex- 
tract the full information. 

The complexity of hadronic r final states re- 
quire sophisticated algorithms based on multi- 
variate technics applied after event reconstruc- 
tion. The next section describes the r reconstruc- 
tion at the D0 experiment. The subsequent sec- 
tion discusses multivariate analysis (MVA) dis- 
criminant, which allows separation of r from jet 
fakes. 



2. TAU RECONSTRUCTION 

In this section, we describe how the r candidate 
is built from elementary reconstructed objects 
within the calorimeter and tracking systems [5]. 
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2.1. Definition of reconstructed objects 

Three objects are potentially used to build re- 
constructed r object: 

1. calorimeter cluster found by a Simple Cone 
Algorithm in a AR < 0.5 cone. 

2. electromagnetic subcluster found using a 
Nearest Neighbour Algorithm with a seed 
in the third layer of the electromagnetic 
calorimeter (i.e. the layer containing the 
finest segmentation). The clustered energy 
is required be greater than 800 MeV. 

3. All the tracks found by the D0 tracking al- 
gorithm in a AR < 0.3 cone around the 
calorimeter cluster are considered. An in- 
variant mass requirement is applied to keep 
only tracks compatible with the r decay. 

In the next section we discuss the use of these 
objects to optimize the r reconstruction accord- 
ing to its hadronic decay. 

2.2. Definition of type candidate 

The D0 experiment bases its r lepton recon- 
struction on specific signatures of the hadronic 
system present in the r final state. Three types of 
candidate are defined and depends on the recon- 
structed objects signature. There are categorized 
as: 

• type 1: one calorimeter cluster and one 
track (optimized for — » n^v T decay); 

• type 2: one calorimeter cluster, one track 
and one calorimeter EM subcluster (opti- 
mized for r ± — » p(— » 7r°7r ± )^ r decay); 

• type 3: one calorimeter cluster, at least 
two tracks (optimized for 3-prong r ± — » 
7r ± 7r T 7r ± ^ r decay). 

2.3. Reconstruction efficiency 

Figure [T] shows the reconstruction algorithm for 
type 2 r. One can see that almost all real r above 
20 GeV are selected by the D0 reconstruction but 
a large fraction of jets (or fake r) remains. In- 
deed, jets from direct strong interaction can easily 
mimic hadronically decaying r lepton and there- 
fore, sophisticated methods is needed to separate 



r from jets. This is the subject of the next sec- 
tion. 
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Figure 1. Reconstruction efficiency for type 2 can- 
didates versus its transverse momentum. Left: 
candidates are genuine r coming from different 
simulated processes. Right: candidates are jets 
from an enriched QCD data sample. 



3. DISCRIMINATION FROM JETS 

3.1. Problematics and strategy 

Hadronic final states of r leptons can be faked 
by jets arising from direct QCD interactions. 
However, jets signatures from r are narrower and 
have lower tracks multiplicity. These properties 
can be used to reject QCD jets and thereby in- 
crease the r identification efficiency. In fact, other 
specificities of hadronic r decays can also be ex- 
ploited and the full information can be combined 
in a Neural Network (NN) [6]. In particular, 
a total of 12 observables (depending on the r 
type) are considered based on a)the isolation in 
the tracking system as well as in the calorimeter 
b) shower shape and composition c) correlations 
between the tracking system and the calorime- 
ter. Figure [2] provides an example of two such 
NN input variables. 

The NN output peaks at for the QCD jets 
and at 1 for the r jets as shown in figure 03 The 
Z — » tt electroweak process is clearly sizeable 
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Figure 2. Example of two discriminating observ- 
ables used in the NN to separate jets from r. Left: 
the observable describes which energy fraction is 
taken by the neutral or the charged component 
of hadronic final state. Right: the observable de- 
scribes if the deposited energy is spread over sev- 
eral calo towers: r tend to have spread deposit to 
jets. 
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Figure 3. The NN output for type 1 and 2 candi- 
dates in fi + r final states. The green histogram 
is the Z — » rr simulation (peaking at 1 as ex- 
pected), the blue histogram is an estimation of 
QCD background and the black dots are the data. 



and data (black dots) are well understood by the 
SM prediction (filled histograms). 

3.2. Identification efficiencies 

Although r reconstruction algorithm is effi- 
cient, the signal over background ratio (S/B) is 
poor. In fact, the order of magnitude of S/B is 
2 at the reconstruction level and reaches almost 
70 after the NN selection. Table [T] summarizes 
efficiencies for each r types after reconstruction 
and subsequent NN identification. 

3.3. Further optimizations 

The optimization strategy is based on how 
the NN treats the information provided by the 
physics. By denoting X = (xi, £2, ■•, x n ) as a 
point in the discriminating variables space, the 
NN output r] nn (X) converges to the true likeli- 
hood function (which is the best classifying func- 
tion): 



rj true (X) 



S(X) 



S(X) + B{X) 



(1) 



Here, S(X) and B(X) are the probability density 
function of the genuine r (signal) and QCD jets 
(background) respectively in the discriminating 
observables space. 



One can find two strategies to optimize the sep- 
aration between jets and r. A first approach is 
to improve rf rue (X) by adding more observables 
based on physical properties of r leptons. A sep- 
arate method would be to minimize \rf rue — r] nn \ 
by performing the NN training in specific phase 
space regions or by increasing the statistics of 
the training samples to describe the differences 
in more details. 

In practice, the first approach was tested 
by trying to include information from the D0 
preshower system [5] which has a better segmen- 
tation than the calorimeter and allows an im- 
proved separation of the 7r° from the 7r ± (for type 
2 decay). After several studies, it has been shown 
that this approach does not significantly help to 
separate r from QCD jets. New physical informa- 
tion can be used such as the long lifetime of the r 
lepton which can lead to displaced tracks from the 
primary vertex similar to the case of b-jets. Fig- 
ure 2] shows a new discriminating variable based 
on the impact parameter of each track for type 3 
r candidates. This new observable increases the r 
identification efficiency of about 10% for the same 
jet rejection. 

A second approach consisting in helping the 
NN convergence was applied by using larger train- 
ing samples and tuning specific NN parameters 
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Table 1 



Fraction (%) of jets and r passing the reconstruction and a NN selection for each r types [7]. 



r type 


1 


2 


3 


all 


jets after reco. only 


2 


12 


38 


52 


r after reco. only 


11 


60 


24 


95 


jets having NN> 0.9 


0.06 


0.24 


0.80 


1.1 


r having NN> 0.9 


7 


44 


16 


67 




(S/B) 

ne 

(S/B) t 



: type all 



Figure 4. Left: signal and background distribu- 
tion of an observable based on the track impact 
parameter (in blue for the QCD jets and in red 
for the real r). Right: a scheme of experimental 
signature of long lived particle. 



such as the numbers of nodes, epochs and mini- 
mization algorithm [6 . Moreover, the difference 
between r and QCD jets can evolve with the can- 
didate transverse energy. In order to help the 
NN to exploit this behaviour, dedicated trainings 
were done for events with high (> 45 GeV) and 
low (< 45 GeV) pr r candidate. 

Overall, the r/jet discrimination is improved 
by 15% at low pr and 40% at high px as shown 
in figure [5] 

4. TAU ENERGY MEASUREMENT 

In order to perform physical measurements in r 
lepton final states, the identification of r leptons 
is necessary but not sufficient. One needs to have 
a good understanding of r kinematic properties. 
This section deals with the r energy measurement 
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Figure 5. Impact of optimizations on hadroni- 
cally decaying r leptons identification versus the 
p T candidate. S/B = N(r)/N(jets) and the 
mention "old" (resp. "new") refers to an event 
selection based on NN without (resp. with) opti- 
mizations. The relative gain is around 20% and 
is better at high pt as expected thanks to the 
specific training. 



at the D0 experiment. 

4.1. Problematics and strategy 

The object energy calibration at collider exper- 
iments is usually performed with the help of well 
known physical processes. At high energies, the 
production and decay of Z bosons allow to per- 
form such calibrations. Two important challenges 
arise for the hadronically decaying r. Firstly, the 
Z — >> rr visible mas£| peak is broad and shifted by 



lr The visible mass observable for Z — >• rr events is defined 
by M^ is = (p Tl + p T2 + Pt) 2 , where ^ is the transverse 



the neutrinos present in the r decay which makes 
this observable less sensitive than for the Z — >• ee 
decay. Secondly, the various r decay modes im- 
plie lower statistics compared with the Z — >• ee 
decay (~ 10 3 events versus ~ 10 5 events). These 
situations require to find another strategy to cal- 
ibrate the visible energy of hadronic r lepton. 

If one considers a type 2 candidate having vis- 
ible energy from (7r°,7r ± ), the true energy can 
be obtained from the measured energy using the 
hadronic calorimeter response (R n ) and the elec- 
tromagnetic calorimeter response (R e ): 



td Tpmeas . z~> rpme 



(2) 



However, the D0 calorimeter is not compen- 
sated (R e R n ) and its segmentation does 
not allow to separate the 7r ± shower from the 
7T° one. Thus, the only measurable energy is 
fimeasj^fimeas ^ j n orc [er to deal with these experi- 
mental constraints, one must use the track energy 
as reference and propagate it to the calorimeter in 
two ways: the absolute correction and the relative 
correction. Both methods are described below. 

4.2. Absolute correction 

The idea of the absolute correction is to use 
the tracker to measure the tt^ energy. In order to 
avoid double counting, the average of the energy 
deposited by charged pion in the calorimeter is 
subtracted: 



E c 



E 



trk 



E cal - {R^(E trk ,ri))-E ] 



■ptrk 



(3) 



Figure [6] shows that method brings a better 
measurement of the r visible energy. 

4.3. Relative correction 

Another strategy to propagate the track energy 
is to correct the simulation event-by-event using 
the E cal /E trk (= E/p) distribution for Z H> rr 
events selected in data: 



(E/p) 



data 



(E/p), 



(4) 



Here, (E/p) is the average value over the events. 
Such a method is relative due to the fact that the 
simulation is adjusted to the data but the true 




missing energy. 



Figure 6. The measured visible energy versus the 
true energy with (right) and without (left) ab- 
solute correction. This study done on simulated 
events show that the absolute correction give a 
better energy measurement. 



energy is not known either in the data or in the 
simulation. However, the R n measurement is not 
needed for this approach. 

5. CONCLUSIONS 

In high energy hadronic collisions, r leptons re- 
quire more sophisticated tools than electrons and 
muons because of their various hadronic decay 
modes and the escaping neutrino. Such r objects 
are of particular interest to probe the Standard 
Model as well as physics beyond the SM includ- 
ing many areas such as the electroweak symmetry 
breaking origin or supersymmetric scenarios. In 
spite of the various experimental challenges, the 
D0 experiment has developed an effective algo- 
rithm to identify hadronically decaying r leptons 
and has performed several physics measurements 
and searches employing such technics. In addi- 
tion, some recent progress has been achieved on 
the r identification improving the potential sen- 
sitivity of the experiment with r leptons. 
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